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This paper is a direct continuation of an earlier work, where we studied Erdos-Renyi random 
graphs perturbed by an interaction Hamiltonian favouring the formation of short cycles. Here, we 
generalize these results. We keep the same interaction Hamiltonian but let it act on general graphs 
with uncorrelated nodes and an arbitrary given degree distribution. It is shown that the results 
obtained for Erdos-Renyi graphs are generic, at the qualitative level. However, scale-free graphs are 
an exception to this general rule and exhibit a singular behaviour, studied thoroughly in this paper, 
both analytically and numerically. 
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I. INTRODUCTION 

The best understood random graphs are those with 
uncorrelated nodes and a local tree structure. It was 
natural to develop the statistical mechanics of random 
networks starting with this highly idealized picture, as 
it is natural to begin explaining the properties of gases 
starting with the ideal ones. However, short loops show 
up in most real networks relatively frequently and it is 
evident that current models should be upgraded to cap- 
ture this common trait (the present state of art in net- 
work research is excellently reviewed in refs. 0,l3il3|)- A 
possible strategy consists in adding to the Hamiltonian 
a term favouring the appearance of specific motifs. This 
strategy was adopted in our recent publication , where 
Erdos-Renyi graphs were perturbed by adding to the ac- 
tion a term proportional to the number of triangles @. 
The purpose of the present paper is to outline a possi- 
ble generalization of the discussion of ref. B| to the case 
where in zeroth order the graphs belong to the statisti- 
cal ensemble of simple (i.e. non-degenerate) graphs with 
uncorrelated nodes and an arbitrary degree distribution 
Pk- 

The partition fmiction of the perturbed model is writ- 
ten in the form 



N 



M 3 = 1 



(1) 



where N and L are the number of nodes and links, respec- 
tively, the sum is over adjacency matrices M, kj is the 
degree of the j-th node satisfying TrlM"^) = kj and 
S{M) is the perturbing Hamiltonian. Setting S{M) = 
one obtains a standard model of networks with uncorre- 
lated nodes and the degree distribution (provided the 



relation '^Pfc/ SfcPfc = 2L/N holds, see for example 
refs 0,0)- In the second line we denote by (• • •) the aver- 
age taken in the unperturbed ensemble and therefore Zo 
is, obviously, the partition function of the unperturbed 
model, in our context an irrelevant overall factor (it can 
be calculated analytically in the large N limit, see eg. |^, 
but we do not need this result here). 

One can, of course, go over from (^| to a grand- 
canonical ensemble, with fluctuating L, multiplying ^ 
by an appropriate L-dependent weight and summing over 
L. In the Erdos-Renyi theory there is a very natural 
and simple recipe for the weight. Hence, in |j] we were 
grand-canonical. Here, we find it more elegant to keep 
the number of links fixed. Notice, that we are interested 
in the large N limit where, for any reasonable choice of 
the weight, L stays anyhow close to its average value as 
a consequence of the constraint L ~ ^ kj and of the 
central limit theorem. 

As in ref. [i| we set for definiteness S{M) = ^Tr{M'^). 
Expanding the exponential in |^ we obtain a perturba- 
tive series representing the partition function. This per- 
turbative representation was thoroughly studied in ref. 
4] and the analytic arguments were completed by nu- 
merical simulations. Let us recall the salient conclusions 
of this study: 

We have introduced a diagrammatic representation 
of possible contributions to the perturbative expansion, 
reminiscent of Feynman diagrams used in field theory. 
Each diagram is a specific subgraph of the full random 
graph. We have shown that the number of contributing 
diagrams grows factorially with the order of the pertur- 
bation theory. 

As is well known, such a factorial growth of the number 
of diagrams signals a breakdown of the perturbation the- 
ory. Indeed, our numerical simulations indicate the pres- 
ence of a transition from a smooth, perturbative regime 
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to a crumpled phase, where almost all nodes form a com- 
plete clique 8] . The transition point Gout {N) scales like 
IniV. The two phases are separated by a "barrier" which 
becomes impenetrable when N ^ oo. 

Remarkably enough it is possible to sum up the lead- 
ing diagrams, i.e. those whose contribution survives in 
the limit TV — > oo. Thus, for example, we have been 
able to derive a closed analytic expression for the av- 
erage number of triangles. It turns out, that at large 
enough A'^ this analytic formula is a very good approx- 
imation in the almost whole region G < Gout{N). Set- 
ting G = GqIuN < Gout (N) one obtains a network 
model with a nontrivially behaving clustering coefhcient 
C cx N^"~^. This clustering coefficient is never constant, 
it falls to zero as ^ oo, but this fall can be made fairly 
slow by a proper choice of Go . 

The results summarized above are the starting point 
of the present paper, which is a direct continuation of 
ref. 3. In Sect. II we extend our diagrammatic rules to 
general uncorrelated graphs with a given degree distribu- 
tion. We also show that the results of ref. continue 
to hold, only slightly modified, in this generalized set-up. 
The so-called scale-free graphs are the only exception to 
the generic behaviour and require a separate discussion, 
presented in Sect. III. Analytic results are confronted to 
numerical simulations in Sect. IV. We conclude in Sect. 
V. 



II. EXTENDING THE DIAGRAMMATIC 
RULES 

The perturbation series is defined as in 4] . One calcu- 
lates the successive terms in the expansion 



([Tr(Af3)]^ 



(2) 



This boils down to the calculation of the expectation 
value of strings Uke Ma^aa^aaaa • • • Ma3,^_^Q3„. A string 
does not vanish iff all the matrix elements it involves are 
equal to unity. The only problem is that the same ele- 
ment, or its transpose, can appear several times in the 
same string and one needs, therefore, to catalogue all 
possible string structures (c/. ref. 0). This is done with 
the help of diagrams: 

Each matrix element Mat, is represented by a line seg- 
ment with endpoints a and h. A string is then represented 
by a collection of n triangles, possibly glued together. In 
order to calculate the perturbation series one has to con- 
sider all possible diagrams. The meaning and the con- 
struction of diagrams is explained at length in ref. 
with the help of explicit examples. We will not repeat 
this discussion here, referring the reader to the original 
paper. We wish only to insist on the salient steps: 

(a) All n-th order terms of the series are, of course, 
proportional to the common factor G"/6"n!. 

(b) Every diagram is a subgraph which has to be em- 
bedded in the full graph. An n-th order diagram has, say. 



V vertices and v < 3n. These v vertices can be identified 
with graph nodes in Nl/{N — -y)! ~ A^^ manners. 

(c) The same diagram topology usually represents a 
number of distinct strings. The calculation of this num- 
ber is the relatively difficult part of the game. But it is 
universal, in the sense that it does not depend on the de- 
gree distribution. The calculations of ref. |i| were done 
for the Erdos-Renyi graph, but hold quite generally. 

(d) Finally, one has to find the expectation value of 
the string corresponding to a given diagram. In Erdos- 
Renyi theory this is simple: if a diagram has £ edges, then 
there are £ independent adjacency matrix elements in the 
string and the expectation value is just p^, where p is the 
control parameter equal to p — a/N when the average 
degree is set to be finite. The case of a general model 
of uncorrelated graphs with a given degree distribution 
requires some extra thought. We will use an argument 
which is not quite original, since it has already been em- 
ployed by other authors in a somewhat different context 
(see, for example, refs. 0,0, 0|) : 

The probability that nodes a, b are connected is 



Prob{a, b) — 



{k)N 



(3) 



where ka{ki,) is the degree of a-th (b-th)node and (fc) — 
J2k ^Pkl TlikPk- Notice that {k)N equals the number of 
directed links. The probability in question is inversely 
proportional to the total number of directed links and is 
proportional to the degrees of the nodes. Eq. © holds 
when the right-hand side is small enough, i.e. at large 
N . The probability that b is in turn connected to, say, c 
is however 



Prob(b, c I a) 



(fcf, - l)kc 
{k)N 



(4) 



because one link emerging from b, the one connecting b 
to a, is already occupied: only fc;, — 1 links are poten- 
tially " active" . Pursuing the argument one derives the 
following rules : 

(a) a factor l/{{k)N) is associated with every edge of 
the diagram, and 

(b) a factor kj/{ka — may. is associated with every 
vertex, say a, of the diagram: here ka is the degree of the 
a-th graph node and ma is the degree of the same node 
regarded as the diagram vertex. 

This was for a particular graph. Averaging, one is led 
to replace 



ka^ika - maY. ^ (fca!/ {ka - ma)!) 



(5) 



in the above rules given above 

We have assumed here that node degrees are uncor- 
related, which strictly speaking is only true in the limit 
N ^ 00. Even in a so-called uncorrelated network some 
"kinematic" correlations appear when one imposes the 
constraint that there are no self and multiple connections 
between nodes and when k"^ /N is not always negligible 
(71 . The last condition is easily satisfied when the degree 
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distribution is defined on a finite support, but may be 
jeopardized in scale-free networks. 

Let us check that for Erdos-Renyi graphs one gets the 
result of ref. 0|: when the degree distribution is Poisso- 
nian the average of a binomial moment is just a power of 
a : 

(fcj/(ft,-ma)!) (6) 
Of course, one has 

Y,ma = 2e (7) 

a 

Furthermore, in Erdos-Renyi theory (fc) = a. Hence the 
diagram with £ vertices gets a factor 

a^'/{aNY=p' (8) 

This is exactly what one has in the grand-canonical 
Erdos-Renyi ensemble and also what one expects in our 
setup in the large limit. 

We are now equipped to calculate the contribution of 
an arbitrary diagram. For example, in the limit G —^ 
the average number of triangles in a graph, the derivative 
of the free energy with respect to G, equals 

- 6 (fc> ) 

because the diagram has three vertices of order two. 
Hence, each vertex contributes a factor (fc(fc — 1)) - re- 
member that vertices are independent - while each link 
contributes a factor 1/ (fc) . The powers of N cancel, as in 
the calculation of ref. Q and the only difference is that 
now instead of a appears a ratio of binomial moments of 
the degree distribution pk ■ 

Higher order diagrams call for other binomial mo- 
ments. It is important to realize that as long as all the 
moments of pk are finite - and therefore iV-independent 
for large enough network size - the hierarchy of diagrams 
in the expansion is the same as in ref. Q. In par- 
ticular, the same diagrams are leading. Summing the 
leading diagrams, i.e. those whose contribution remains 
finite when N —i- oo, one gets 




which generalizes eq. (36) of ref. Higher order mo- 
ments (cumulants) of the T-distribution are obtained by 
differentiating the right-hand side of with respect to 
G. One finds that the distribution is Poissonian. 

The situation changes when the moments of pk are not 
necessarily finite. Without much loss of generality we will 
limit our attention to the so-called scale-free graphs, i.e. 
those where pk falls at large k like a power: pk oc . 
We also assume that (fc) is finite: (3 > 2. 




FIG. 1: Additions of a triangle to a diagram, leading to the 
increase of the number of diagram links by three (a), two (b) 
and one (c). 



III. SCALE-FREE GRAPHS 

When Pk ^ a/k^ at k ^ 1, the moments of order 
larger than /3 — 2 diverge. Actually, at large but finite N, 
the degree distribution of a non-degenerate uncorrelated 
graph is cut at 

kma.'xN^ 7 = min(l/2,l/(/3-l)) (11) 

(see ref. for a derivation of this result). Hence, higher 
order moments of the degree distribution increase like 
some powers of N and, consequently, the hierarchy of 
the expansion is modified compared to ref. It 
turns out that three cases require a separate discussion: 



A. Case /3 > 3 

We will prove that the leading diagrams are the same 
as those considered in ref. 4] , i. e. those where the num- 
ber of triangle edges is equal to the number of triangle 
vertices. One can easily see that in these diagrams the 
triangles can overlap, but otherwise do not touch. 

The proof is by induction. First we observe that the 
diagram of order n = 1 (one triangle) is leading and its 
contribution tends to a constant as — > cx3. Higher 
order diagrams can be constructed by adding successive 
triangles. Suppose that at the n-th order the leading 
diagrams are those with the number of links equal to 
the number of triangles. What happens when one adds 
the {n -\- l)-st triangle ? If it is put on top of an existing 
triangle or if it is isolated, the new diagram scales with N 
the same way as its ancestor. In general, the new triangle 
is not isolated and its addition increases the number of 
links. Three possibilities, illustrated in Fig. 1, can arise. 
Because of (lllll increasing the degree of a diagram vertex 
by Am links produces at most an extra divergent factor 

Nf^ (12) 

When the number of links is increased by three (Fig. la) 
Am = 2 in three vertices. The new links yield a factor 



4 




FIG. 2: Leading diagram for P — 3. 



N ^ each. Thus the new global factor is at most 



N- 



9-3g 



(13) 



Similarly, for Figs. lb,c one finds that the potentially 
dangerous new factor is at most 



6 — 2f3 3-/3 

N -s-i and TV^-i 



(14) 



respectively. For /3 > 3 these factors do not diverge. We 
conclude that the leading diagrams are those considered 
in ref. Q. The average number of triangles is finite in 
the limit N —> oo and is given by eq. 110(1 . 



As explained in ref. P|, calculating the average num- 
ber of triangles at large N it is sufficient to consider con- 
nected diagrams only. At each order of perturbation the- 
ory the leading connected diagram is the one with the 
largest number of vertices of order 2, i.e. the diagram 
shown in Fig. 2. In the n-th order we have n triangles 
and therefore the number of vertices of order 2 is 2n. 
The central vertex has order 2n. Each triangle in Fig. 2 
can be glued to the central vertex in three manners. The 
contribution of the n-th order diagram to the free energy 
is therefore (for n > 1): 



1 



-(fc(fc-i)r 



6"7i! 2n- 2 (fc)3 



In^" N 



2(n-l)7i! (8(fc)3)" 
where we have used the the asymptotic result 

(fc(fc- 1)) - ^IniV 
In this limit eq. © reads 



(19) 



(20) 



B. Case /? = 3 

As already recalled, at finite N the degree distribu- 
tion is cut and the cut-off scales like ^/N. Of course, the 
cut-off is not sharp, but in calculations it is convenient 
to replace it by a sharp one kmax — c\/N. The deter- 
mination of the constant c is not obvious. We make a 
choice which yields more or less correct moments of the 
cut degree distribution: at large enough N 



(fc^)«;^lniV 



and 



(15) 



(16) 



(the sharp cut-off upsets the normalization, but the effect 
is negligible for large N) . The constant c can be roughly 
estimated using eq. (flHjl : 



In^iV 



(21) 



Summing all the connected diagrams one gets the free 
energy. Differentiating the latter with respect to G one 
finds the average number of triangles, which can be writ- 
ten in the form 



where Tq — {T)g=o and 



63 



(22) 
(23) 



It is evident that the increase of (T) with Grj* cx Gin N 
is nearly exponential. 



N 



(17) 



Let us calculate the leading iV-dependence of an arbitrary 
diagram with £ edges and v vertices: one has a global 
factor (we omit logs) 



(18) 



since each link yields N~'^, the v vertices yield and 
the binomial moments associated with v vertices yield the 
sum in the exponent. However, this sum can be calcu- 
lated: —2) = 2£— 2v. Hence, the total exponent 
is zero, only logs remain. 



C. Case 2 < /3 < 3 

For /3 < 3 all the moments, apart from the first, diverge 
like a power of . As shown in 7] the cut-off of the degree 
distribution scales again like ^/N. Therefore an arbitrary 
diagram's contribution is proportional to the factor 



N 



^-i+iJ2^irn,-fi+l) ^ ^§(3-/3) 



(24) 



Hence the leading diagrams are those with a maxi- 
mal number of vertices. One can calculate the num- 
ber of triangles (T) in the leading-power approxima- 
tion. Again, it is sufficient to consider connected dia- 
grams only. Assuming here, for simplicity of writing, 
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TABLE I: The average number of triangles at G = : {T}g^o 
from Monte Carlo simulation (left) and estimated from degree 
distributions using eq. @ (right). 



AT 
i V 


o or; 

p — Z.O 
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p — o 


O A 

P — 4 


1024 


26.36(9) 35.49 


11.01(3) 13.15 


4.91(2) 5.31 


2048 


42.26(18) 54.91 


15.10(7) 17.58 


5.80(2) 6.14 


4096 


67.05(36) 84.16 


20.26(12) 22.93 


6.58(3) 6.85 


8192 


106.2(9) 129.5 


27.22(28) 30.28 


7.33(3) 7.56 


16384 


167.0(1.2) 198.9 


35.54(33) 38.80 


8.01(9) 8.19 


TV ^ oo 


48 


i In-^ N 

5 
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that (fc™) = c™7V2(™+i f3) g^^^ using the combinatorial 
coefficients calculated in ref. we get: 

It is obvious that the result is very sensitive to the value 
of c. This would also be true if we were using a more 
realistic ansatz, viz. the analogue of (|16|l . A much more 
reliable prediction is the approximate scaling law: 



(T) -ro/(GTff) 
In the large N limit eq. reads 

To^(T)G=oociVi(3-/3) 



(26) 



(27) 



and f{x) : f(0) = 1 is a positive, monotonically increas- 
ing function. Although we have calculated only the first 
few terms in the expansion H25(l , an educated guess is that 
the rise of the right-hand side is again nearly exponential. 



IV. NUMERICAL RESULTS 



in ref. ji| we use the algorithm of ref. , which has the 
advantage of generating not only nondegenerate graphs 
but also thermal fluctuations. The latter point is im- 
portant because our main goal is to check the stability 
of the smooth, perturbative phase. We will also com- 
pare the numerical Monte Carlo data to the predictions 
of the preceding section, obtained by summing leading 
diagrams. The agreement will be only semi-quantitative, 
because of large finite-size corrections. 

Actually, three types of effects occur at finite N: First, 
eqs. Q-© are only approximate, especially when the 
degree distribution has a fat tail, as already explained. 
Second, nonleading diagrams are negligible, at fixed G, 
only asymptotically and it turns out that asymptopia 
is hard to reach. Furthermore, the nonleading contri- 
bution blows up as one approaches the transition point 
G = Gout, because the transition is an intrinsically non- 
perturbative phenomenon. These two effects are in a 
sense a nuisance for us, since they just obscure the pic- 
ture. Third, there are manifestations of the fat tail 
Pk ~ in the node degree distribution that persist 
at any N and are, in fact, increasingly important as N 
grows up. They are responsible for the singular behaviour 
of scale-free networks and are therefore of much physical 
interest. 

As in our earlier publications we assume, for definite- 
ness, that the degree distribution has the form (c/. ref. 

M) 

Similar results were obtained with other choices of pk- 

In Table 1 is given the average number of triangles at 
G — 0, viz. {T)a=o, measured directly in Monte Carlo 
simulation (the left column) and estimated using eq. © 
(the right column) from the degree distribution generated 
in the same simulation for /3 = 2.5, 3 and 4, respectively. 
The agreement improves as N increases and worsens, as 
expected, with decreasing P and growing fat tail. Notice, 
that {T)g=q increases with N also for (3 = 4, although it 
is expected to be finite in the limit N oo: the second 
moment of the degree distribution is still rising signifi- 
cantly in the explored range of N. 
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In this section we continue the study of the preceding 
one but having recourse to Monte Carlo simulations. As 



FIG. 3: (r)/r„ versus G for /3 = 4 and iV = 2"' (x), 2" (+), 
2^2 (V), 2^^ (A), 2" (O). 




FIG. 4: ((T) - To)/Ti 



versus GT^ 



(x),2" (+), 2^^ (V), 2" (A), 2^ 



for (5 

(O). 



3 and iV 



In Fig. 3 we plot {T) /Tq versus G for /3 = 4 and 
N ranging from 2^° to 2^'^. The asymptotic expectation 
exp(G') is also drawn. The data scale and follow the curve 
exp(G') at small G. The characteristic fan shows up as 
one moves toward the transition point G — Gout- No- 
tice, that largest N data tend to be closer to the curve, 
as expected. At this point we are unable to tell how 
Gout behaves. We recall that Gout oc InA^ when in ze- 
roth order the graphs are of Erdos-Renyi type. This is 
expected to be the generic behaviour when the degree 
distribution has effectively a finite support. In this re- 
spect the status of the /3 = 4 case is uncertain. There 
is no evidence in the data for an increase of Gout with 
TV. The apparent constancy of Gout could be a finite-size 
effect, however. On the other hand, it is very plausible 
that hubs, i.e. nodes with largest degree, behave as seeds 
of Strauss cliques, preventing the growth of the barrier 
separating the smooth phase from the crumpled one. In 
order to settle this question one would have to simulate 
enormous networks, beyond reach with present means. 

In Fig. 4 we plot ((T) - To)/t| versus Gr|, as sug- 
gested by eq. for /3 = 3 and N ranging again from 
2^" to 2^^. The data scale reasonably well, especially at 

2 

low GTq . It appears that, very roughly 



Gout ~ ^ cx In ^ 

3 

^0 



(29) 



The constant c can be estimated from eq. (|17|) . using 
the observed (i.e. cut) degree distributions. A reliable 
estimate is obtained from low order moments and using 
rather large data; one finds c ranging from 0.9 to 1.2. 
Setting c = 0.9 one gets from eq. H22II the curve shown 
in Fig. 4. 

2 

In Fig. 5 we plot (T) /Tq versus GT^ , as suggested by 
eq. H25|l . for /3 — 2.5 and A^ ranging from 2^" to 2^*^. One 

observes the expected scaling at small values of GT^ , but 
at larger values finite-size effects become gradually more 



FIG. 5: 

2" (+), 



(T)/To versus GT^f for p = 2.5 and N = 2^° (x) 
2^2 (V), 2" (A), 2^^ (O). 




1000 



2.5. Dashed line is for G = 0, solid line for G = 0.08 (just 
before the transition to the crumpled phase). The almost 
straight line corresponds to the asymptotic shape of Pk, i.e. 
to Pk. 



and more important. It appears that, very roughly 



2.3 

Gout ~ — oc N 

T 3 
^0 



(30) 



In Fig. 6 we plot the degree distribution at A^ = 2^^ 
and P = 2.5, for G = and 0.08 ( just before the tran- 
sition to the crumpled phase), to show that the intro- 
duction of the interaction Hamiltonian does not spoil the 
scale free property. In Fig. 7 we show the variation with 
k of the clustering parameter 



Gfc 



fc(fc- 1) 



(31) 



where Tk is the average number of triangles touching a 
vertex with degree k. One observes that the decrease of 
the clustering parameter is much slower than observed in 
Internet, for example, where Ck oc k~^-'^^ [l7| . 



V. SUMMARY AND CONCLUSION 

This paper is a direct continuation of ref. 4], where 
we have studied perturbed random Erdos-Renyi graphs. 
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FIG. 7: Clustering coefficient versus node degree at A'' = 
2^'*,/? = 2.5 and G — 0.08 (just before the transition to the 
crumpled phase). Notice, that the scales are linear. 

In the present paper we show that the results of ref.jj 
continue to hold when an ensemble of (almost) arbitrary 
uncorrelated graphs is perturbed by the same interaction, 
favouring the formation of triangles. There is, however, a 
notable exception to this generic behaviour: the so-called 
scale-free graphs behave differently, especially when the 
degree distribution has a diverging variance. 

At finite N the smooth phase exists only when the 
interaction coupling G is smaller than some threshold 
G < Gout- Gcnerically Gout scales like InA^, but for 
scale-free networks with 2 < /3 < 3 it scales like N^~^ 
(modulo logs). Hence, the support of the smooth regime 
does not expand but shrinks to zero in the thermody- 
namic limit. There is nothing dramatic in such a be- 
haviour, which is also encountered in more conventional 
matrix models. It just means that G is not the physi- 
cal coupling and that the latter is rather GN^^^, again 
modulo logs. When the physical coupling is used the 
smooth phase lives in a finite coupling interval and the 
average number of triangles is enhanced compared to the 
unperturbed expectation. However the stability of this 
phase is an open problem. In ref. P| we found that 
nonperturbative phenomena are manifestly negligible al- 
most in the whole smooth phase: numerical data were 
remarkably close to perturbative predictions. Here, we 
cannot claim the same, partly because we are unable to 
disentangle finite-size and nonperturbative effects. 

At finite G scale-free networks are unstable for N large 
enough. The physical significance of this singular be- 
haviour is not fully clear yet. In our simulations the ther- 
mal motion consists of network rewirings. Rewiring is an 
ergodic move: every two states can be transformed one 
into another by making a finite number of rewirings. In 
particular, any other algorithmic move could be regarded 



as made up of rewirings. However, with a different algo- 
rithm the thermalization time of the system would in 
general be different, in particular it could explode with 
increasing N. Thus, it is not excluded, although does 
not seem very likely, that the instability is an algorithm 
artifact. 

Very many natural networks are scale-free, with the 
exponent f] below 3. At least some of them seem fairly 
stable. In some cases there are selection rules constrain- 
ing the rewirings. But this is not the most interesting 
possibility. How do the natural networks compare to the 
graphs of our model? We see one significant difference: in 
our graphs the clustering coefficient is weakly correlated 
with node degrees, while in natural networks it tends to 
decrease like some power of the latter. The behaviour 
of our graphs is easy to understand as follows: The sys- 
tem forms triangles at random and therefore tends to 
attach many triangles to hubs, which apparently become 
seeds of Strauss cliques. Natural networks seem to avoid 
this disease by suppressing the formation of triangles at 
hubs. It is plausible that a specific hierarchical organi- 
zation screens natural networks from the instability. In 
such a scenario a triangle generating term with a finite 
coupling G could after all be present in the Hamiltonian. 
For the moment this is just a speculation. 

Our paper is not a phcnomcnological one. We are 
not yet at the stage of constructing a model to be com- 
pared to the data. We focus on the theoretical problem 
of the stability of networks with respect to motif gener- 
ating terms in the Hamiltonian. However, the paper is 
not quite devoid of phenomenological implications: our 
method allows us not only to calculate averages of phys- 
ical quantities characterizing an individual network, but 
also fluctuations of those quantities in the ensemble, giv- 
ing us an insight into the problem of typicality of net- 
works. As far as wc know the magnitude of fluctuations 
of motifs has never been estimated for graphs with an 
arbitrary given degree distribution. 
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